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Abs tra ct 

This paper discusses the trends in longitudinal 
dynamic aeroelastic stability of a generic X-Wing aircraft 
model with design parameter variadons. X-Wing rotor 
blade sweep angle, rado of blade mass to total vehicle 
mass, blade structural stiffness cross-coupling and vehicle 
center-of-gravity locadon were parameters considered. The 
typical instability encountered is body-freedom flutter 
involving a low frequency interacdon of the first elasuc 
mode and the aircraft short period mode. Parametric cases 
with the lowest stauc margin consistendy demonstrated 
the highest flutter dynamic pressures. As mass rado was 
increased, the flutter boundary decreased. The decrease was 
emphasized as center-of-gravity location was moved 
forward. As sweep angle varied, it was observed that the 
resulting increase in forward-swept blade bending 
amplitude reladve to aft blade bending amplitude in the 
first elasuc mode had a stabilizing effect on the flutter 
boundary. Finally, small amounts of suffness cross- 
coupling in the aft blades increased flutter dynamic 
pressure. 

.Ipaa riu ct i oa 

The X-Wing aircraft is a unique vehicle conflguradon 
combining the verucal take-off advantages of a helicopter 
and the high-speed forward flight capability of a fixed 
wing aircraft. The aircraft concept uulizes a four biaded, 
bearingless rotor system capable of operauon in three 
flight phases. Verucal take-off and hover with rotadng 
blades comprise the first flight phase. The second phase 
is a conversion mode in which the blade rotadon slows to 
a stop. In the third phase, the stopped rotor blades act as 
fixed wings for high speed forward flight [1]. Two of the 
blades are swept forward and two are swept aft 
symmetrically. From this conflguradon, the aircraft has 
become known as the "X-Wing.” 

In each of these phases of flight, circuladon control 
provides a mechanism for augmenting and controlling lift 
[2]. Essentially, circuladon control is achieved by 
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blowing a sheet of air through spanwise slots over the 
trailing edges of Coanda airfoils. Because these quasi- 
ellipucai airfoils have rounded trailing edges, the flow 
remains attached and the airfoils' stagnadon point is 
relocated. Thus, lift and stability can be modulated by 
controlling the velocity and spanwise distribudon of 
blowing. 

It has been demonstrated in the past that fixed-wing 
aircraft with forward-swept wings exhibit substandal 
coupling between elasuc and rigid-body vibration modes. 
This symmetric, low frequency flutter condition is 
commonly referred to as 'body-freedom' flutter and is a 
coupling of wing bending and rigid-body pitch and plunge 
motions [3], Because of the similarities between this 
conflguradon and the X-Wing in a fixed rotor mode, it is 
likely that such an instability will be characteristic of the 
vehicle. 

Several studies have focused on the aeroelastic 
stability of this unique configuration. For instance, in 
1981, Gimmestad [4] conducted a study of an X-Wing 
configuration which included rigid-body freedoms as well 
as circuladon control blowing along the blades. The 
results indicate that aft biade motion damps unstable 
forward biade motion through body freedoms until high 
velocities when the vehicle loses stability in pitch and a 
whole vehicle divergence results. 

In 1987, Gilbert and Silva [5] showed analytically 
that with increasing airspeed aeroelastic deformations of an 
X-Wing configuration caused a forward shift in the aircraft 
aerodynamic center (a.c.) location. The result, again, is a 
loss of static margin or vehicle divergence. Compared to 
the previous study, however, the divergence occurred at 
lower airspeeds which is due to lower blade natural 
frequencies and higher pitch inertia. It was also concluded 
that no antisymmetric divergence or flutter modes exist. 

More recently, Haas [6] has investigated an X-Wing 
configuration with circulation control blowing but 
without rigid-body freedoms. At high angles of attack, a 
single degree of freedom flutter involving first bending 
elastic mode occurs. The instability is due to the airfoil 
section’s negative lift curve slope at high angles. It is 
emphasized that the flutter is not dynamic stall, as there is 
no flow separation. At reduced angles of attack, classical 
bending-torsion flutter is observed. With no circulation 
control blowing, the highly rigid vehicle experiences 
classical bending-torsion flutter and static divergence at 
very high airspeeds. 

The intent of this study is to determine the dynamic 
aeroelastic behavior of the generic X-Wing aircraft 



configuration with body freedoms. As a basis for the 
study, a generalized X-Wing aeroelastic model is developed 
for the vehicle in a stopped rotor mode. Composite-beam 
finite elements were used to model the X-Wing blade 
structure and Doublet Lattice lifting surface theory [7] is 
used to calculate the unsteady aerodynamics. Variations in 
both structural and aerodynamic parameters are made to 
determine trends in the flutter behavior. The analysis uses 
only symmetric vibration modes and does not include 
circulation control blowing. 

Aeroelastic Vfodel Development 
Equations of Vfotion 

The equations of motion for a free- flying aeroelastic 
vehicle are. in terms of vehicle rigid body and elastic 
vibration modes [8], 

H Ml s" + [C ]s + [K ] + 

q[Q(s)]]{nCs)} = 0 0) 

where q is the dynamic pressure, [QC s>] is the matrix of 
generalized aerodynamic forces, {*H( s ) } is the vector of 
generalized coordinates, and s is the Laplace transform 

variable, [m'1 , [c'\ and [k'J are the generalized mass, 
damping and. stiffness matrices, respectively. 

Structur al Mo del 

To model the generic X-Wing configuration, a 
structural half-model was developed using beam finite 
elements. The model was general in the sense that 
parametric variations could easily be made in the sweep 
angle, mass and stiffness properties. The beam element 
was incorporated into EAL (Engineering Analysis 
Language) for free vibration analysis of the model [9]. 

The finite element used in this study was developed to 
model beams demonstrating bending-torsion stiffness 
cross-coupling and is described by elemental mass and 
stiffness matrices which are fully defined in reference [10]. 
Both elemental matrices involve a nondimensional 
parameter, V, which describes the stiffness cross- 
coupling. It is defined by Weisshaar [ 11] as 


y=± V if /El • GJ (2) 

where El and GJ are the bending and torsional stiffnesses 

of the beam, respectively, and R denotes stiffness cross- 
coupling between bending and torsional deformations. 
Limits on V are derived from the energy requirement that 
a stiffness matrix must be positive semi-definite 
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or equivalently 

EICJ-Kf >0. 0 

which implies 

- i < y < i 

A positive value of ¥ indicates that upward 
(positive) beam bending induces a leading edge up twist or 
'washiri. A negative value implies negative twist with 
positive beam bending or 'washout'. ¥=0 represents the 
absence of stiffness cross-coupling. Physically, ¥ is 
determined by the composite ply orientation within a 
given composite laminate. 

The actual X-Wing structural half-model takes 
advantage of the vehicle's symmetry while in the fixed 
rotor mode. It consists of two untapered blades, one with 
forward sweep and one with aft sweep, each modeled by 
ten finite elements. The specific nodal layout and some 
parameter definitions are illustrated in Figure 1. 

Vehicle center-of-gravity (c.g.) location, also shown 
in Figure 1, is determined by the two mass components of 
the X-Wing half-model. One mass component is the set 
of two blades. The other is an attached mass simulating 
the fuselage and is positioned at varying locations along 
the vehicle's longitudinal axis. The element used to link 
fuselage and blade motion is massless and rigid. 

Some typical natural frequencies and modeshapes 
computed in this study are shown in Figure 2. 

A ero dyn amic Mo del 

The X-Wing half-model is further developed by 
defining aerodynamic lifting surfaces. Combining these 
pianform geometries with the free vibration modeshapes 
and frequencies, an unsteady aerodynamic analysis can be 
conducted to determine generalized aerodynamic forces 
(GAFs) acting on the vehicle at various reduced 
frequencies. 

Two assumptions have been made in the aerodynamic 
calculations. First, because symmetric structural 
vibration modeshapes are used for analysis, the 
aerodynamics are also specified as symmetric about the x-z 
plane. Second, circulation control blowing has not been 
included. Although it was shown by Haas [6] that the 
strength and spanwise distribution of blowing directly 
affect the static and dynamic aeroelastic vehicle responses, 
it is assumed in the present study that there is no 
circulation control along the blades. 

Unsteady GAFs in this investigation are computed 
using the Doublet Lattice Method as available in ISAC 
(Interaction of Structures, Aerodynamics and Controls) 
[12]. They are tabulated as functions of Mach number and 

reduced frequency, k, defined as k= where ca is the 

oscilladon frequency, b is one-half the blade chord length 
and V is the velocity. GAFs are computed for harmonic 
motion and are extended to arbitrary motion using a 
rational function approximation method described below. 
In these computations, Mach=0.3 and sea level conditions 
are assumed. A representative Doublet Lattice panel 
layout is shown in Figure 3. 


Tabic I Nominal parameter values 


The rational function approximation to the tabulated 
GAFs is made in the form [8] 
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where L is the number of aerodynamic lags and 0 . are 
the aerodynamic lag coefficients. [a J, 


and the 


B. 


's are real and are computed using a least 


squares fit to the tabulated GAFs for the case of harmonic 
motion, s=j co. The lags, 0 . , are arbitrarily chosen to be 

within the reduced frequency range under investigation and 
provide a good approximation for s=j to. The 
approximations are constrained to be an exact fit to the 
tabulated GAFs at k=0 in order to define as accurately as 
possible forces generated by rigid body modes and the 
steady state aerodynamics. In this study, there are four 
aerodynamic lag coefficients, 0 = 0.1, 0.15, 0.25, 0.4 and 
ten values of reduced frequency in the range of k=0.0 to 
k=1.0. 


Slate Sp ace F ormu lation 

With the aerodynamics in functional form and the 
structural matrices defined, the equation of motion, Eqn. 
1, can be rewritten in state space form for stability 
analysis as shown in the Appendix [8]. 


A 

45° 

sweep angle 

m r 

0.25 

blade-to-vehicle mass ratio 

V 

0.0/0.0 

stiffness cross-coupling 



in forward/' aft blades 

X 

-0.15 

center-of-gravity locadon 


The IS AC stability analysis for each case involves 
determining the eigenvalues of the associated system 
matrix [F], Eqn. A-6 in the Appendix, over a range of 
velocities. A complex eigenvalue or root is stable when 
the real part is less than zero. A positive real part implies 
unstable, i.e. positively damped, oscillatory motion of the 
associated rigid-body or elastic mode. A real root is stable 
if iL is less than zero and unstable if it is greater than zero. 

NQmina l C on fi gurati on 

To illustrate typical trends in X-Wing stability, the 
root locus of the nominal forward c.g. case ( x =-0.15) is 
presented in Figure 4. Velocity variation is in approx- 
imately 30 mph intervals between 30 and 600 mph. 
Since only the first two elastic modes and rigid body 
modes show any tendency towards instability, higher 
modes have not been shown. First and second elastic 
mode roots repel one another; the first mode frequency 
decreasing while the second mode frequency increases. 
Coupling of the first elastic mode with the vehicle short 
period mode creates a body-freedom flutter condition. 

Maa-BatiQjy.anatioQ 


Aerp ei as tic Stabili ty . Anaiasa 
A na lysi s . Me thod 

Having developed a generalized X-Wing model, 
parametric variations may now be made and the vehicle 
stability studied. A parametric variation is the changing 
of some vehicle dimension, stiffness or mass property 
while fixing the other system parameters to some nominal 
values. Each specific combination of parametric and 
nominal values is referred to as a 'configuration'. In this 
study, center-of-gravity, sweep angle, blade- to- vehicle 
mass ratio and stiffness cross-coupling are chosen to be 
the system parameters. The nominal value of each 
parameter is listed in Table I. 


Mass ratio, m r , is defined for the half- model as the 
ratio of total blade mass over total vehicle mass. 


where mf is the attached fuselage mass and mb is the mass 
of two blades. The nominal value of m r is 0.25. Pitch 

inertia as well as total vehicle mass are held constant 
during m f variations. Configurations are examined at 

m r =0.15, 0.25, 0.35, 0.45, and 0.55, in each case at two 

forward c.g. locations, x =-0.15 and x =-0.3. 

Velocity root loci are, for most of the cases, similar 
to those shown in Figure 4. First and second elastic 
mode roots move away from one another as velocity 
increases; interaction between the short period and first 
elasdc roots results in the latter roots moving to the right- 
half plane. The exception is for the case with m r =0.15 

and x =-0.15. Contrary to the typical trend, as velocity 
increases, the short period roots are driven to the right-half 
plane while the first elasdc mode roots are forced to the 



left. The instability mode, however, is still body-freedom 
flutter. 

Several interesting results were obtained from the 
stability analysis. These results are summarized in Figure 
5 which presents the flutter dynamic pressure as a function 
of mass ratio. Figure 5 shows that at m^O.15, flutter 
boundaries for x =-0.15 and x=-0.3 cases intersect. 
Further investigation shows that these two cases are 
essentially the same configuration. With small mass 
ratios, the fuselage is massive enough to create a 'near- 
clamped' fuselage condition so that c.g. location is no 
'longer a factor. Furthermore, since the structural 
characteristics of the blades are the same for a given m r at 
any c.g. location, the first two elastic mode natural 
frequencies are nearly identical. 

As mass ratio increases, the effect of c.g. location is 
to separate the flutter boundaries. Cases with the more 
forward c.g. locations lead to the lower flutter dynamic 
pressures. In addition. Figure 5 shows that an overall 
decrease in flutter dynamic pressure, qf, occurs as mass 
ratio increases. This happens because blade first mode 
natural frequencies decrease by as much as 50% when 
is increased from 0.15 to 0.55, as shown in Figure 6, 
allowing short period and first elastic modes to interact at 
increasingly lower dynamic pressures. The decrease occurs 
because blade mass per unit length, which is inversely 
proportional to the square of natural frequency, increases 
with m r 

Sweep Angle Variation 

As shown in Figure 1, blade sweep angle is defined as 
one-half the relative angle between forward- and aft-swept 
X-Wing rotor blades. The nominal value is A= 45°. 
Configurations examined are A = 15°, 30°, 45°, 60°, and 

75°, each with c.g. locations x = -0.15 and -0.3. 

Velocity root loci for the various sweep 
configurations are similar to those shown in Figure 4. 
The first elastic mode and short period mode roots interact 
to create a body- freedom flutter condition. However, as 
A increases, a transition occurs. The A =15° cases 
have short period mode roots moving to the right half 
plane at flutter, whereas, the A =45° cases have first 
elastic mode roots moving to the right half plane. This 
transition is shown in Figure 7, a root locus for the 
configuration described by A =30° and x=-0.3. There is 
a velocity at which first elastic mode roots and short 
period mode roots become indistinguishable. It is unclear 
which roots pass into the right half plane. A similar 
transition also occurs in the x=-0.15 cases but at a 

sweep angle not investigated (15°< A <30°). 

Structural vibration modeshapes for representative 

cases A = 15° and 45°, both with x=-0.l5, are compared 
in Figure S. In the First elastic mode of the A =45° case, 
both blades move together with the aft blade having a 

larger relative amplitude. As A decreases below 45°, a 

"modeswitch'’ occurs and, in the extreme case of A =15°, 
it is seen that forward and aft blades move opposite to one 


another although the aft blade still has a larger relative 
amplitude. Thus, at some transition sweep angle, denoted 
as A tr, the first elastic mode involves, primarily, bending 
of the aft blade and no relative forward-blade bending. A tr 
increases with forward motion of the c.g.. For x =-0.15 

cases, A =* 30° while for x =-0.3 cases, A =* 45°. 
tr tr 

Results of the dynamic stability analysis are shown 

in Figure 9. Several interesting observations were made. 

In both x=-0.15 and -0.3 cases, minimums in flutter 

boundaries are present. The sweep angle at which these 

minimums occur are denoted by A m . In the x=-0.15 

case, A ^ * 38°and in the x =-0.3 case, A _ » 45°. A 
m m 

qualitative investigation of the static margin variation 

with A has shown that these minimums are associated 

with maximums in static margin [13]. This is consistent 

with previous results obtained. In the m r variations, it 

was shown that increasing static margin led to decreasing 

flutter dynamic pressure. In addition, as the c.g. location 

is moved aft, the maximum static margin occurs at a 

higher sweep angle and, thus, the stability curve is shifted 

to the right. The stability boundary minimums are 

associated, too, with the modeshape transition. For the 

most unstable case, the first elastic modeshape involves 

very little or no forward-blade bending amplitude relative 

to aft-blade bending amplitude. The a.c. is in its most aft 

position and static margin is maximum or near 

maximum. As forward blade bead Lag amplitude increase s 

or decreases with variations in A , the a.c. is shifted 

forward and the static margin is reduced. Thus, forward- 

blade loads counteract alt-blade loads and tend to stabilize 

the vehicle. 

Aeroelastic Tailoring 

The final objective of this study was to find a means 
of delaying or eliminating body-freedom flutter of the X- 
Wing. For this purpose, aeroelastic tailoring through the 
use of stiffness cross-coupling was investigated. This 
coupling was incorporated independendy into the forward- 
and aft-swept blades of a nominal parametric configuration 
and in varying amounts: ¥=-0.5, -0.25, 0.25, and 0.5. 

Velocity root loci for the tailored X-Wing 
configurations were found to be very similar to the typical 
root locus shown earlier in Figure 4. All cases experience 
body-freedom flutter as a result of short period and first 
mode interaction. 

The vehicle flutter dynamic pressure as a function of 
blade stiffness cross-coupling is shown in Figure 10. In 
this figure, one flutter boundary is associated with cases 
in which ¥ changes only on the forward blades and 
remains zero in the aft blades. The other curve represents 
a boundary for cases with ¥ changes in the aft blade 
while forward-blade cross-coupling is zero. It is seen that 
stiffness tailoring of the forward-swept blades has a small 
effect on the flutter stability of the parametric case 
evaluated. The result is not surprising because the 
forward-swept blades have a small bending amplitude 
relative to the aft-blade bending amplitude in the first 
elastic mode of this configuration. Thus, changes in 
aerodynamic loading due to coupling-induced washrn or 
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washout of the forward blades are small when compared to 
the aft-blade loads. 

With negative cross-coupling, 'washout', in the aft- 
swept blades, the natural tendency of the blades to 
washout during bending is magnified. As a result of the 
decreased aerodynamic loads, the static margin is decreased. 
It has been shown previously that a decrease in static 
margin increases flutter dynamic pressure. The addition of 
positive cross-coupling, washin, in the aft blades causes 
the flutter dynamic pressure to decrease below the 
uncoupled flutter dynamic pressure. The natural tendency 
' of the blade to washout has been dominated by stiffness- 
induced washin. As a result of the increased blade loading, 
static margin is increased and flutter dynamic pressure 
declines. 

There exists a point of diminishing returns, where 
further increases V lead to decreases in qf. This may be 
explained by considering the decrease in first mode natural 
frequencies as the magnitude of ¥ increases. The 
frequencies change slightly because both the mass matrix 
elements and the bending stiffness are functions of ¥. 
As the magnitude of ¥ increases, variations in the 
mathematical model increasingly affect the system 
characteristics. In this case, the lower frequencies allow 
for interaction of short period and first clastic mode roots 
at lower dynamic pressures or qp 

P ncn .MQm^ni_Effg,£is 

An analysis involving vehicle stability as a function 
of c.g. location yields some interesting insight into 
vehicle stability trends. The locations considered were 
x=-0.3, -0.15, -0.075, 0.0, and +0.15. Figure 11 
summarizes the results. The upper curve defines a 
boundary above which body-freedom flutter occurs. For 
x<-0.15, flutter dynamic pressure is slightly less than 
that of the nominal case. It increases dramatically as the 
c.g. is moved aft, which is effectively decreasing the static 
margin. The lower curve is a restabilization boundary of 
the vehicle divergence condition which is consistently 
observed in the cases where the initial a.c. location is 
forward of the vehicle c.g.. 

As illustrated in the root locus of an initially unstable 
case. Figure 12, the short period mode root lies along the 
positive real axis and moves outward with increases in 
velocity. The root eventually reaches a limit and begins 
to the left along the real axis. Restabilization of the root 
is followed at higher velocities by the body-freedom flutter 
condition described previously. 

Insight is gained by considering the first elastic 
modeshape as a function of c.g. location. Figure 13. In 
the forward c.g. locations, aft-swept blade bending 
amplitude is much larger than forward blade bending 
amplitude. At x =0.0, blades bend together with the same 
amplitude because the vehicle is longitudinally 
symmetric. As c.g. is moved aft, forward blade bending 
becomes the dominant bending amplitude. Herein lies the 
possibility of vehicle divergence for the statically unstable 
cases. At low velocities, forward swept blades begin to 
develop larger airloads due to aeroelasuc washin. Because 
forward blade bending amplitude is significantly larger 


than aft blade bending amplitude in these cases, there is 
insufficient nosedown pitch moment generated by the aft 
blades to counteract noseup moments due to the forward 
blade airloads. Thus, the X-Wing experiences a vehicle 
divergence instability characterized by an aperiodic upward 
pitch. As velocity increases, forward-swept blade inertial 
loads [15] and aft blade airloads increase, theoretically 
restabilizing the vehicle. One implication of this is that 
stability augmentation of the divergent motion may only 
be required at low velocities. 

Similarly, in the statically stable cases, forward-swept 
blade loads serve to counteract unstable aft blade loads. At 
the most forward c.g. location, the modeshapes involve 
large aft blade bending amplitude relative to forward blade 
bending amplitude. These configurations demonstrate 
the lowest flutter dynamic pressures. As c.g. moves aft, 
there is an increase in the relative forward blade bending 
amplitudes as well as a decrease in the longitudinal static 
margin, both of which cause reductions in the vehicle 
pitch moment. In addition, the body-freedom flutter 
dynamic pressures improve. Thus, the same decrease in 
vehicle pitch moment which encourages vehicle 
divergence discourages body-freedom flutter. 

This statement is supported by aeroelastic tailoring 
results. When aft blades experience cross-coupling 
induced washout, an overall decrease in pitch moment 
occurs. Figure 10 emphasizes the associated increase in 
flutter dynamic pressure. 

It is noted that the tendency to have an increasing 
body-freedom flutter dynamic pressure with a decreasing 
longitudinal static stability margin has been shown in two 
previous studies. One study involved flutter of the X-29 
[3] and the other involved flutter of a tailless sailplane 
[15]. As another point of interest, the addition of a tail 
would probably increase nose down pitching moment and 
therefore would shift the dynamic divergence curve to the 
right as shown in Figure 11. The flutter boundary, 
however, would drop because increases in pitch moment 
have been shown to be destabilizing. 

Conclusions 

The X-Wing’s typical aeroelastic instability while in 
the fixed-rotor mode is body-freedom flutter, a low 
frequency interaction between the first elastic mode and the 
short period mode. Over c.g. variations, parametric cases 
with the negative static margins demonstrated the highest 
flutter dynamic pressures. As the ratio of X-Wing rotor- 
to-fuseiage mass increased, the flutter boundary decreased. 
This trend is due to a drop in first elastic mode natural 
frequency with increasing mass ratio. The decrease in 
flutter dynamic pressures was emphasized as c.g. location 
moved forward. An increase in forward-swept blade 
bending amplitude relative to aft-swept blade bending 
amplitude in the first elastic mode due to sweep angle 
variations had a stabilizing effect on the flutter boundaries 
at two different c.g. locations. Finally, negative stiffness 
cross-coupling or washout when incorporated into the aft 
blades caused the flutter dynamic pressure to increase 
above the uncoupled' flutter dynamic pressure. Positive 
cross-coupling in the aft blade or washin was 
destabilizing. In the forward-swept blades, stiffness cross- 





coupling had little effect on stability for the parametric 
case evaluated. Finally, it was concluded that in most 
cases forward blade loads and aft blade loads tended to 
counteract each other in the presence of body freedoms. 

Appendix 

Substituting the aerodynamic approximation, Eqn.3, 
into Eqn. I yields 
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where 

[m]= [Ml +q[A 2 j(y) 


In the time domain, the state space form is 


z= [F]z 


(A-7) 


[t]=[Cl +q[Aji 
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First order state equations for the aerodynamic 
approximation terms can be obtained as 
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Fig. 2 Modeshapes and frequencies for nominal 
case. 
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node locations. 
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Fig. 5 Flutter dynamic pressure as a function of 

m and x . 
r 



Fig. S Variation of First elastic mode as A increases 
for configurations with x = - 0. 15 . 



Fig. 6 Variation of natural frequency with x and m r 
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in sweep angle. 
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Fig. 7 Root locus of A = 30° . x = - 0. 3 

configuration illustrating interaction of 
short penod and first elastic modes. 



Fig. 10 Variation of flutter dynamic pressure with 
V as it is applied independently to the 
forward-swept blade tVf) and to the aft- 
swept blade (vytj). 
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Fig. 11 Stability boundaries as function of x, 
including the effect of a tail. 
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Fig. 12 Root locus illustrating re stabilization of 
the short period mode roots along the real 
axis. 


V 







i i 


*0.3 -0.15 0.0 0.15 0.3 

x 


Fig. 13 Transition of First elastic mode vertical 
displacement as a function of x . 
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